Effective Potential and Thermodynamics for a Coupled Two-Field Bose Gas Model 
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We study the thermodynamics of a two-species homogeneous and dilute Bose gas that is self- 
■ - - interacting and quadratically coupled to each other. We make use of field theoretical functional 

integral techniques and evaluate the one-loop finite temperature effective potential for this system 

■ considering the resummation of the leading order temperature dependent as well as infrared con- 
tributions. The symmetry breaking pattern associated to the model is then studied by considering 
different values of self and inter-species couplings. We pay special attention to the eventual appear- 

■ ance of re-entrant phases and/or shifts in the observed critical temperatures as compared to the 
^ ' mono atomic (one- field Bose) case. 

PACS numbers: 03.75.Hh, 05.30.Jp, ll.lO.Wx, 98.80.Cq 

Q . I. INTRODUCTION 

-J i 

Due to a larger region of parameters, multi-field theories may exhibit a much richer phase structure in comparison 
to single field theories. Phenomena that are completely absent in single field theories may arise in some regions 
^2 ' of the space parameter related to multi-field theories. From a qualitative point of view one expects that, at finite 
H ] temperatures and/or densities, multi-field theories will display phase transition patterns which cannot occur in the 
' ■ single field case (see for instance Refs. 0,0 and references therein). Quantitatively, on the other hand, the actual 
y ' value of critical quantities may be different in the two cases. For example, consider a scalar multi-field theory composed 
Q , by two different fields, $ and which self interact via quartic interactions (e.g., 17$$^ and 5xX*)- The two fields may 
^ ' also interact quadratically with each other, e.g., with an interaction term g^^x^- general, these type of interactions 
lead to an 0(iV$) x 0{N^) invariant theory, where Ng, and represent the number of components of each field. In 
^ the relativistic case one may think of <i> and x as representing an extended Higgs sector of the standard model, e.g. the 
^ Kibble-Higgs sector of a SU{5) grand-unified theory, with iV$ = 90 and = 24 In the non- relativistic case one 
. can associate $ and x to two different species of bosonic atoms in an homogeneous dilute Bose gas with = = 2. 
' To assure stability the numerical values of the couplings need to observe certain constraints whose general form is 
5* 5x ^ ^5^' where the self-couplings satisfy (7$ > 0, > 0, whereas g can be either positive or negative, b is some 
'^O ' positive real number whose value depends on the way the interactions are normalized in the Lagrangian density. At 
the same time, stability (boundness of the potential energy) only requires that the coupling constant be positive in the 
one field case. It is then easy to understand what was said above about the numerical values of critical quantities, like 
the critical temperature, Tc, which in many cases depends on the masses and couplings, causing T™"'*' jmono^ Also, 
"j^ ' the way the transitions occur can be highly influenced by the presence of the crossed interaction term proportional to 
I g, as first noted by Weinberg This fact arises as a consequence of 5 assuming either positive or negative values and 
I , still keep the theory bounded from below. In his work, Weinberg used perturbation theory to analyze a relativistic 
' 0(iV$) X 0{Ny^) model at finite temperature. He found that, for g < 0, one may find regions of the parameter 
Ch , space where unexpected phenomena arise. In particular, it was found that a symmetry that was broken at T = 
^ ■ could remain broken at arbitrarily high temperatures, in what was called symmetry nonrestoration (SNR). Also, a 
' symmetry that was unbroken at T = could be broken at some finite Tc in a manifestation of inverse symmetry 
^ 1 breaking (ISB). It is worth recalling that these two phenomena never show up in the mono-field 0{N) theories where 

one always reaches the symmetric phase at some finite temperature. 
^ As phenomena happening at finite temperatures, one could argue that SNR and ISB are just artifacts of perturbation 

theory, which although used in Ref. "T], is well known to be inadequate to treat high temperature field theories 5, 6]. 
However, powerful non-perturbative techniques P,Q)I3 that include resummations of leading order and infrared terms 
and thermal effects on the couplings, have confirmed the possibility of these phenomena showing up in the relativistic 
case. The nonrelativistic case, where the appearance of such phenomena is highly counter intuitive, has been treated 
very recently by some of the present authors 0, 0, There, it was shown that SNR and ISB cannot manifest 
themselves when thermal effects on the couplings are taken into account. A broken symmetry at T = is always 
restored at some finite Tc while an unbroken symmetry (at T = 0) can be broken at a finite T^^^ only to be restored 
at a higher critical temperature characterizing a re-entrant phase that is typical of many condensed matter systems, 
like the Rochelle salt, spin glasses, compounds known as the manganites, liquid crystals and many others materials, 
as recently reviewed in |l]| . In the work carried out in Ref. one body terms proportional to (and k^X^) have 

been introduced to drive symmetry breaking. Then, apart from the masses and couplings one has the parameters k 
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which may, for example, represent external fields that are temperature independent at the tree level. 

The aim of the present work is also to analyze a multi-field nonrelativistic theory but this time we want to make 
contact with the case that is relevant for dilute homogenous weakly interacting Bose gases. Therefore, we will choose 
a U{1) X U{1) version (or equivalently, 0{2) x 0{2)) taking the one body parameters as representing the chemical 
potentials. As we shall see this choice makes the finite temperature treatment very complex and here we will use 
functional techniques to evaluate the effective potential (or free energy) . This will allow us to investigate the possibility 
of SNR/ISB occurring in coupled homogeneous Bose gases. At the same time we will be in position to check if the 
coupling of two different species of Bose gases shifts the critical temperatures in relation to the case where this 
coupling does not exist. Due to the non-perturbative nature of such evaluation we will consider a resummation of the 
leading order temperature dependent and infrared contributions to the effective potential to one-loop. As we shall 
see, the present investigation excludes the possibility of re-entrant phases when the nonrelativistic models studied 
in Ref. Q contain the chemical potentials needed to represent Bose gases. Moreover, the values of the critical 
temperatures for the two different species seem to be insensitive to the existence of the (new) cross coupling. As in 
the mono atomic case, their values only depend on the density (and mass) of each specie coinciding with the ideal 
gas situation. Although the conclusion about the non existence of re-entrant phases is very plausible the statement 
about the critical temperature should be taken with care. This is because here we are only resumming one-loop (or 
direct contributions) which do not contribute at the critical point due to the Hugenholtz-Pines theorem. On the other 
hand, our non-perturbative calculation can reveal the differences in the thermal behavior of uncoupled and coupled 
Bose gases at least in the region < T < T^. The point T — can be fully exploited only by using non-perturbative 
techniques which go beyond one-loop [l^ . 

The experimental realization of the Bose-Einstein condensation (BEG) in dilute atomic gases has greatly stimulated 
an enormous number of theoretical studies in this field (for recent reviews on the theory and experiments, see for 
instance Ref. 0|). Most of this interest comes from the fact that in these experiments a great deal of control can be 
achieved in almost every parameter of the system. Thus, experiments in dilute atomic Bose gases provide a perfect 
ground to test numerous models and field theory methods applied to these models, as for example finite temperature 
quantum field theory methods as commonly used to study phase transition properties of relativistic models. This 
makes BEG one of the most attractive and promising systems in which one can use models and can test schemes and 
approximations that could also prove useful in very different environments such as in the early universe, heavy-ion 
collisions, etc, where the use of one or multi-field models may have importance in their understanding. Therefore, the 
extension to BEG systems of the analysis performed in previous works in the phase structure of multi-field models is 
particularly interesting given the possibility of using experiments in dilute atomic Bose gases as an analog system to 
model and test these finite temperature quantum field theory systems. 

This work is organized as follows. In Sec. II we present the model representing two, self-interacting, coupled Bose 
gases. In Sec. Ill we evaluate the finite temperature effective potential to one-loop order. Section IV, which is divided 
in two subsections, presents the evaluation of thermodynamic quantities as well as self-energies. In the first subsection 
we consider the one-field case whereas the two-field case is considered in the next subsection. In Sec. V we consider 
the high temperature (symmetry restored) phases. Our conclusions and final remarks are presented in Sec. VI. 



II. THE TWO SELF-INTERACTING COUPLED MODEL FOR BOSE GASES 



Let us now consider the case of two coupled Bose gases in the presence of self-interactions. The model we consider 
is similar to the ones used in other theoretical studies of homogeneous dilute coupled Bose gases , that consists of 
a hard core sphere gas model described by nonrelativistic interacting (complex) scalar fields, with an overall repulsive 
potential. This system can be described by the following global U^{1) x [/$(!) invariant Lagrangian model for 
two-species nonrelativistic complex scalar field $ and 4*, with self-couplings and and inter-species coupling g, 



£($*,$, I-*,*) = (z5t + ^V2 + -^($*$)2 
\ 2to$ / 2 

+ ^* (idt + ^— ) * + /I**** - ^(***)2 - g{^*<^>){^^'*^^f) , (2.1) 

\ 2mq, J 2 

where the associated chemical potentials are represented by fii {i = 'i> or <I>) while represent the masses. For 
the hard core sphere self-interactions we take the phenomenological coupling constants as being the ones normally 
used in the absence of cross interactions and which are valid in the dilute gas approximation |l3l |. In terms of the 
corresponding s-wave scattering lengths, a^, they can be written as gi = Anai/rrii, while the cross-coupling is chosen 
as g = 47ra*_$/m*^$ where m^_$ = m'},m$ / (jn^ + m$) represents a reduced mass. 
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The use of analog condensed matter systems to study multi-field theory models like (|2.1() , can be envisaged by the 
use of a system composed by a mixture of coupled atomic gases, like the ones originally produced in Ref. [13 in 
which one has the same chemical element in two different hyperfine states and that may be treated as "effectively 
distinguishable" , or just consider the mixing of two different mono-atomic Bose gases. 

The equivalent finite temperature Euclidean (r = it) spacetime action to Eq. H2.1|l is given by 



3„ 



SEif3)= dT d'x 







9 V2 „^^,^7*,^,^^2 



'i>*(^-^-/i$)'i>+Y($*$)'+.9(***)($*<f) 



(2.2) 



where T — l/j3 (we are considering throughout this paper all quantities in natural unities, where ?i = fc^ = 1). Let 
us initially consider '5 and $ as pure condensates, in the absence of fluctuations (pure c- numbers). By varying the 
action Eq. (|2.2|) with respect to ^ and $ we obtain the equations of motion (in Euclidean time), analogous to the 
Gross-Pitaevskii equation in the case of a monoatomic gas, given by 

dr^ = -77^ - + (5*|*|2 + g|$|2) , 

2m* ^ ' 

9.$ = -1^^^'^ - + (5*1*1' + 5l*l') $ • (2.3) 

The minimization of the potential part of the action Eq. (|2.2|l leads to the relations involving the condensates for ^ 
and $, c = I^P and p$.c = respectively, and the chemical potentials: 

M* = 5$P$,c + 5P*,c ■ (2.4) 

In addition, for Eq. I|2.4|l to represent a local minimum of the potential, we still need to make the requirement that 
it be bounded from below, which requires the coupling constants to satisfy 

3*5$ - 5^ > , (2.5) 

with repulsive self-interactions, 17* > 0, 5$ > 0. In the experimental situation of a coupled binary system of Bose 
atoms, the condition H2.5|l is required for the mixture of condensates to be stable. In the case where Eq. (|2.5(l is 
violated, de-mixing of the condensates happens, so that the mixture of condensates will tend to separate spatially, as 
actually observed in the experiment of the second reference in Ref. 15]. 

Equation (|2.4(l also gives the condensate densities in terms of the chemical potentials: 

P*,c 2~ ' (2.6) 

5*5* - 5 

and 

P<E.,c = — ■ (2.7) 

5*5* - 5 

III. THE ONE-LOOP FINITE TEMPERATURE EFFECTIVE POTENTIAL 

Let us now consider the effect of fluctuations. Here we choose to study the system at equilibrium and will include 
fluctuations in the system by means of the field theoretical method of the effective potential. We will follow closely 
the functional integration derivation used by the authors in Ref. to the usual one field self-interacting Bose gas 
(for other field theoretic methods applied to BEG problems, see also Ref. 0,0] and for a review Ref. 0|). As 
usual in the computation of the effective potential, we start by decomposing the fields $ and ^' in Eq. H2.1|l in terms 
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of (constant) background fields (which, without loss of generality, can be taken as real fields, c- numbers) 4>q and ipo, 
respectively, and (q-numbcr) fluctuations $ and which in terms of real components, become 



'i' = (?Ao + V^i + #2) ■ 



(3.1) 
(3.2) 



Considering the leading order in the fluctuations (which is equivalent to keep terms up to order h in the effective 
potential, or one- loop order), when substituting Eqs. I|3.1|l and (|3.2I) in Eq. (|2.1() we only need to keep the quadratic 
terms in the fluctuation flelds for the computation of the one- loop potential for the background fields, c/jq and -00 • We 
then obtain the functional partition function to one-loop order 



= J D^iDiP2D4>iDct>2 exp (-52) , 



(3.3) 



where the functional integral is restricted over fields satisfying the periodic boundary conditions (the Kubo-Martin- 
Schwinger condition) ^'(x, t) = ^'(x, r + /3) and <f>(x, r) — $(x, t + (i). 5*2 is the Euclidean action to quadratic order 
in the fluctuation flelds, 



S2 



dr I drx 



'^'^ , 5* , 4 A^* , 2 , 5* , 4 , .9 .2 / 2 , 1 n> 



(3.4) 



where we have defined the vector x — (^ij '^2, "011 "02) and M is the matrix operator for the quadratic terms in the 
fiuctuations, 



idr 



M 









iVf. _ ^ J- ,/,2 I 9^2 








(3.5) 

The partial time derivative in (|3.5() is over Euclidean time: dr = d/dr, r = it. As usual, the effective potential is 
defined from the functional partition function by 



149(00,00) = -— lnZ[ 



(3.6) 



where V is the volume. By performing the functional integration in the quadratic fluctuations x, the one-loop effective 
potential Vcs{4'0i4'^) obtained from Eq. (|3.4(l is then given by 



Kff(0o>o) = ~^4>l + - ^00^ + f + + ilndetM , 

where the last term on the RHS of Eq. H3.7|) comes from the functional integral over the components of x, 



(3.7) 



ilndetM^-^ J 13011)02 £'V'ii502 exp dr J d^x (^^x ■ M ■ x 

Expressing Eqs. H3.8|l and H3.5|l in the space-time momentum Fourier transform form, we obtain 



(3.8) 



■IndetM 



1 1 

2P 



E 



n— — 00 



d3q 



(27r)= 



9' 



w*(q) 



5* ,2 , .9 ,2U 
^0o + 20oJ| 



(3.9) 



5 



with 



and 



(3.10) 



where LOi{q}, z = <i>, ^' is given by 



'^<&(q) + — 00 ' 2 



tj$(q) 



(3.11) 



w»(q) 



q 

2m, 



(3.12) 



and in Eq. I|3.9|l represents the Matsubara frequencies, u;„ — 27rn//3, n = 0,±1,±2, • • •. 

Considering the classical condensates densities in the absence of fluctuations, Eqs. (|2.6(l and (|2.7|l . noticing that 
P'i'.c = s-nd /5$^c = </'o/2 and substituting them in Eqs. H3.10|l and (|3.11|) we recover the Bogoliubov dispersion 

relations for the gases in the broken phase 



and 



[w^(q) + 35*p*,c + ff/o-i>,c] [w*(q) + ff*P*,c + 



/ q2 



2m5( \2m4( 



(3.13) 



[u}^{q) + 'ig<i,p^,c + 5P*,c] [w*(q) + 54>P*,c + 5P*,c] 



2m$ \2m$ 



(3.14) 



which are consequences of the breaking of the two continuous symmetries of the model and Goldstone's theorem. We 
can then recognize in Eqs. (|3.13|l and H3.14(l a Higgs and Goldstone modes like terms, that we denote respectively by 
Hi and Gi, given by 

(3.15) 
(3.16) 
(3.17) 
(3.18) 

In terms of Eqs. (|3.15|) - (|3.18|l . the contribution of fluctuations to the classical potential in (|3.7|l . AV, becomes 

n ^ ' 

d3q 



iJ*(q, ( 


/>O,'0o) 


= w*(q) - 


3.9* ,2 


+ 2^0 


Gv[,(q,( 


/•o, V'o) 


= t^*(q) - 


.9* ,2 , 




i7$(q, ( 


/•o, V'o) 


= t^$(q) H 






G$(q, ( 


/>O,'0o) 


= w$(q) H 




1^0. 



1 1 

2^ 

1 1 
2^ 



E 



(27r)3 

d3 



ln[K + A2) K + S^)] 



(3.19) 
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where the terms A and B in Eq. (|3.19() are given by 



A' 



and 



2 H^G^ + H^G^ 1 



B' = 



(i7*G* - H^G^y + Ag^^''^<i>^G^G^ 



1/2 



1/2 



(3.20) 



(3.21) 



The one-loop correction to the classical potential when expressed in the form of the last term in the RHS of Eq. (|3.19fl 
is a suitable form that allows to easily perform the sum over the Matsubara frequencies by using the formula, 



1 ;^ 1 1 



(3.22) 



Defining the quantity v{lu), 



,2 , , ,2^ 



(3.23) 



we have 



which, from Eq. 13.22|l . becomes 



and so 



dvim) _ 1 



- y 



2uj 



n— — oo 



dv[ijj) 
du! 



= 1 



(3.24) 



(3.25) 



v[ijj) = cj + — In (l — e ^'^) + terms independent of 

p 

Neglecting the constant terms, we then obtain for Eq. H3.19|l the result 



UJ 



(3.26) 



(3.27) 



IV. THERMODYNAMIC QUANTITIES AND SELF-ENERGIES 

Given the effective potential, from Eqs. (|3.7() and (|3.27|l . 



— *» 



5* ,4 



■9j,2 ;2 



(4.1) 



with the functions A = A(q, 0oi'0o) and B = B {q, (j)Q , tpo) , given by Eqs. H3.20|l and (|3.21|l . respectively, we can 
compute all relevant thermodynamical functions. In particular, we have that the pressure is defined as the negative 
of the effective potential computed at its minima (which is the thermodynamical free energy of the system), 
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P = P{T, Ai*) = -Vofr(T, 00, V'o) 



(4.2) 



where (f)„i and tpm are the values of 0o and that extremizes (corresponding to a minimum of) the effective potential, 



dV,s{T, 00, ^o) 



<t<a=<t>m,4'0 = 4'r, 



0, 



ayeff(T,0o,V'o) 



aVo 



(4.3) 



For the tree-level potential, and (j)„i are given by Eqs. (|2.6|l and H2.7|l . respectively (with ^ = V'm/^ and p$_c = 
0^/2). As we will see below, interactions will change these tree-level expressions and, consequently, modifications 
will have to be implemented in the effective potential so as to preserve the Goldstone's theorem in the presence of 
interactions and finite temperature effects. From the pressure the total number of particles follows as 



P* = 7^ , P* = ^ • (4-4) 

The corrections to the tree-level (zero temperature) densities change Eqs. (I2.t)|l and H2.7|) according to Eq. H4.3|) and 
the obtained expression for the effective potential. Let us see this with some more details starting with the uncoupled 
gas case. 



A. The One-Field Case 



For 5 0, Eqs. and (jX^ become 



3.9* 



cj$(q) + 



5* 



(4.5) 



and 



B'{g = 0)= (co^id) + (c.*(q) + , (4.6) 

and the contributions from the $ and 5* fields to the effective potential, Eq. H4.f|l . decouple. Lets consider then, e.g., 
the $ field contribution. In this case 




The zero temperature contribution in Eq. (|4.7() is divergent and require proper renormalization. It is easier to do 
it by performing the momentum integral in d = 3 — e dimensions and the resulting integral is found to be finite in 
dimensional regularization (when taking e at the end). This is so since the divergence in Eq. I|4.7|) is a power-law 
one and, therefore, in dimensional regularization, the regularized integral results to be finite. We then obtain the 
minimum of the effective potential as 



5* 



260$ (q) + ^ 



{2n) 



[I + 2n$(q)] 



(4.8) 
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where 



n$(q) 



exp 



(4.9) 



is the Bose-Einstein distribution for the single self-interacting field and, in Eq. H4.8|l . fJ.<i>/g<t> is the tree- level condensate 
density (which also follows from Eq. H2.7I) for g = 0). By demanding that the spectrum for the single self-interacting 
gas at finite temperature is still gapless (Goldstone's theorem), we can define the condensate density like for instance 
(/)^/2 — fl^/gi^, where p,<^ denotes an effective chemical potential (we here will be using an analogous definition as 
taken by the authors of Ref. in their study of the effective potential for a single self-interacting Bose field) . This 
is expected, since including the finite temperature contributions implies that the original chemical potential must be 
changed accordingly^. In terms of these new definitions Eq. (|4.8|) becomes 



5$ 



2-3- 



(27r)3 



\j 2m4. ( 



2m^ 



2/i$ 



exp 



2m4 



2/j,$ 



(4.10) 



Replacing fi^ by the effective chemical potential in the fluctuation terms of the effective potential corresponds to 
replacing the tree-level field propagators by the self-energy dressed ones. In this case, the equivalent of the matrix 
M, Eq. H3.5|l . for the one-field case, in the basis {(j)i,(t>2) and in momentum space, becomes 



2 



-UJ„ 



^4 

2 ^0 ■ 



01 ,<P2 

9i^2 I ■ 



(4.11) 



■^02,02 



One-loop diagrams that contribute to the self-energies are shown in Fig. ^ 



D 



(a) 



JJ 




(b) 



D 




FIG. 1: Typical one-loop Feynman diagrams contributing to the field self-energy. External lines stand for (j>i and the internal 
propagators Djk are defined by the inverse of the matrix of quadratic fluctuations M. 



Since Y,^.^^- for i ^ j are linear in [j„ (see, e.g. Ref. 0|), they are identically zero, 5^01,02 = 
Eq. (|4.11|) . the equivalent of the Higgs and Goldstone mode terms, Eqs. (|3.17|) and (|3.18|) becomes 



— 0. From 



G<]j 



2m$ 

^2 



3.g$ 
2 

5$ ,2 



00 + 



2m$ ■ 2 ^^"^^ 

Therefore, Goldstone's theorem is preserved if we define the effective potential p,<s> as 



(4.12) 
(4.13) 



The same effect happens in thermal field theory, where instead of a thermodynamic chemical potential, we would then now talk about 
a constant mass term for the field. However, finite temperature contributions entering via the self- energies change this mass such as 
to make it temperature dependent. This is equivalent to change the original mass by a "dressed" one, where self-energy corrections are 
taken into account in the definitions of the field propagators. 
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— /i$ — S02^02 , (4-14) 

and the minimum of the effective potential, 0,„, is related with /!$ by 0^„/2 — Note that at the critical 

temperature of Bose-Einstein condensation, /x$ =0(4>m = 0) and Eq. (|4.14|) becomes /i$ = E^j.^j, which is just the 
standard form for the Hugenholtz-Pines theorem |19l |. 

The introduction of the self-energies can be implemented at the Lagrangian density level right at the 

beginning in a self-consistent way that avoids any possible overcounting of diagrams as would be caused by the naive 
change of the matrix of quadratic fluctuations given by Eq. I|4.11|l . It can be checked that, without any further 
modifications, the change produced by Eq. (|4.11() starts overcounting diagrams at the two-loop level (like two-bubble 
vacuum and self-energy diagrams in perturbation theory) . This overcounting can be eliminated completely by writing 
the original Lagrangian density, in terms of the field components (jji and (/)2 like: 



i 

where the self-energy terms are at the same time added and subtracted in the original action. While the added 
self-energies dress the field propagators, the subtracted terms are treated as additional interaction terms and they 
here act in the sense of subtracting the extra contributions coming from the dressing of the Higgs and Goldstone 
modes, Eqs. (I4.12|) and H4.13|l . The procedure in H4.15|l is common in various other instances of studies involving 
resummation of quantum and temperature correction terms in quantum field theory |2Clj | . 

From the self-energies, one can check that for temperatures T ^ which is the regime of temperatures we are 
interested in, the exchange diagrams (like the two-vertex one-loop diagrams of the form of Fig. lb) are subleading 
compared to the tadpole diagrams. Fig. la. (see also Ref. 16]). In this case, S^j^^^j^ ~ ^02:02 they are given, at 
the minimum of the effective potential, by 



2o^ 

2m^ 



{2n) 



1 + 



2m^ 



exp 



q" 

2m^ 



2m<] 



- 1 



(4.16) 



Comparing Eq. (|4.16() with H4.10|l . we get the result given previously for the condensate density at high temperatures, 
0^/2 = fl^/g^. 

In terms of Eq. (|4.11|) . the dressed effective potential therefore becomes just like in Eq. (|4.7|) . but with 
exchanged by /!$ in the correction terms for the tree-level potential. With this self-energy "improved" effective 
potential computed at the minimum (f>Q — (/).,„, we obtain the pressure, 



P(T,Ai*) = 



25* 
(2^ 



XT) 



1 




(4.17) 



where we have neglected the corrections coming from the new interaction terms given by the change of the Lagrangian 
density as in Eq. (|4.15|) . since these terms are higher order than the ones considered in the one- loop level. From Eq. 
(I4.17|l . we determine the total density, p$ as^ 



^ Note that here and also for the two-field case discussed below we are discarding contributions coming from the derivatives of the self- 
energy, both with relation to the fields and chemical potential, since these terms result to be higher order than the one-loop order being 
considered. For instance dT,/d(j>Q is already of order 0{g^) and same order corrections at the two-loop order should be considered for 
consistency as well. 
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5$ 



(2^)3 



or, also using Eq. H4.1U|I to express in term of 



2^ (2^ + ^Ai* 



exp 



/3 



2^ ( 2^ + 2/^4 



(4.18) 



5* 



rf^q 2^ + 



(27r)3 



2m.^ \ 2m.^ 



2/i,j 



1 + 



exp 



(4.19) 



The T = term in the RHS of Eq. (|4.19fl gives the quantum depletion of the condensate, while the finite temperature 
term gives the thermal depletion, as usual. The T = momentum integral term can be easily computed using 
dimensional regularization. Converting the momentum integral to arbitrary d = 3 — e dimensions and taking e ^ 0, 
we obtain the result, 



1 f 2^+ A* (4.20) 



27 



2m.^ \ 2m 



Writing Eq. H4.19|l in terms of the temperature dependent condensate density, p<s>^c{T) = 0^(T)/2 with = f;$(/)^/2, 
taking it at the critical point, T = T^, and having that fl^(T = T^) = (or ipmiT = T^) — 0, since the condensate 
density vanishes at Tc), we immediately obtain from Eq. H4.19|l that 



J (27r)3 _ ^ LC(3/2)J 

where C(3/2) ~ 2.612. The result Ea. l4.21() is the standard one for the critical temperature for an homogeneous Bose 
gas^. 

It is instructive to see, from Eq. (|4.19|) . how the condensate density p<^,c{T) changes with temperature. Writing 
the coupling constant 5$ in terms of the s-wave scattering length 47ra$ /m$ and defining the dimensionless 

quantities p<j)^c — P<s>,c/P'S>, = T/T^^c, where T^^c is given by Eq. (|4.21|l . we obtain for Eq. H4.19|l . the implicit 
equation for p^^dT), 



37ri/2 



1/2 -3/2 



^1/2^(3/2) 



t; 



3/2 



dx X 



x^ + 2C(3/2)2/3n^/^p^^,/f^ 
x2+4C(3/2)2/3n^/3p^_jf*^'^' 



:[a;2+4C(3/2)2/3„y3-^^/j^^j 



(4.22) 



where we have made the change of integration variable in Eq. 1)4. 19|) . = / (2m^T) and also used the quantum 
depletion result Eq. H4.20|l . In Eq. H4.22|l we have defined the diluteness parameter 71$ = p$a| The integration 
in X in (|4.22|l can easily be performed numerically producing the standard result shown in Fig. 12 for different values 
of n$. 



^ Note that at this level of approximation that we are considering, Tc is the same as that for the ideal Bose gas. Corrections due to the 
self-interactions are only acces sible thro ugh no n-perturbative methods (beyond one-loop) requiring at least second order corrections in 
the self-energy, see, e.g., Refs. |12| and 11911211 for recent reviews. 
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FIG. 2: Results for the dimensionless p$,c as a function of r$, for n# = 10 (solid) ,n$ = 10 (dashed) and n,j, = 10 
(dotted). 



B. The Two-Field Case 



Let us now return to the two-field case model. The field propagators in the non-vanishing self-energy contributions 
are found by the inverse of the dressed matrix of quadratic terms, M, which generalizes Eq. (|4.11(l of the one-field 
case and Eq. I|3.5|l for the free (inverse of) propagator terms in the (0i, 02, "011 ^"2) basis. In momentum space M is 
represented by the matrix 



M$,*(w„,q) 



g(j>oipo + Erf 

V 








g4>o^o + ^01,^1 




Gil 






- E, 



(4.23) 



■l/'2,V'2 



where the functions Hi, Gi, ? = are given by Eqs. (|3.15|l - (|3.18ll . respectively. In obtaining the self-energies 

is easy to show that all cross- like self-energies vanish identically except by S0i ,0i, which is given by exchange- like 
diagrams (like in Fig. lb). The introduction of the self-energy terms in H4.23|l can again be implemented self- 
consistently already at the Lagrangian level by an analogous procedure as the one shown in Eq. I|4.15|l for the 
one-field case. 

In complete analogy with the one-field case we now have, from Eq. (|4.23() . that in order to preserve Goldstone's 
theorem in both and 4* field directions in the broken (condensed) phase, the inclusion of fluctuations must change 
the tree-level chemical potentials to 



= /i* 



(4.24) 
(4.25) 



As in the one-field case, we restrict our main analysis to the high temperature regime, for which T ^ /Ij. In this 
regime it can again be checked that the exchange diagrams are subleading compared to the tadpole diagrams. In this 
case E0i,^i can be neglected compared to the tree level term in Eq. H4.23|l and S^i,^ 
They are given, at the minima of the effective potential {(prmipm), by 



and 



(2^)3 



dA „ 



dB 
Wo 



<l>0 ^ <l>rr 



(4.26) 



(27r)3 



dA 



dB 



4>Q = 4-77 



(4.27) 
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with 



(4.28) 



where A and B, which follows from the determinant of Eq. H4.23|l are the analogous of Eqs. (|3.20l) and (|3.21() with 
/i$ and /Zip in those expressions replaced by the effective chemical potentials /x$ and /x^, given by Eqs. (|4.24|) and 
(|4.25|l . respectively. The minima of the effective potential, (j)m and V'mi are determined from the effective potential 
Vcs{T, 00, V'o), Eq. H4.1|l . with A ^ A and B ^ B. Note that in this case, just like in the one-field case, the condensate 
densities in the presence of fluctuations, p^i^dT) = '0m/2 and p^.c{T) = (p^/'^j are given by the same equations H2.6() 
and 1)2.7(1 . but with chemical potentials given by the effective ones. 



p*,c(r) 



g^g^ - g'^ 



(4.29) 



and 



- g 

Let us now compute the total densities, and which are obtained from Eq. 



(4.30) 

with the pressure given by 



1 



2 (5*5* - 5^) 



(27r)= 



(A + B) 



In 1 



(l-e-'^^) +ln(l-e-'^^) 



1/JQ = 1/J7, 

00 = 'f^n 



(4.31) 



If we now express A as A = A[V'm(M*, A**, 2^), 0m(M*, M*, 2^), M*, A**] and B as B = 
B['4'm{p-<s>, P'S',T),4>mip<s>, P'^,T), p<i,, pqi] it foUows, from Eq. (|4.4II . that we can write for p^ and pq, the fol- 
lowing expressions: 



Pt = 



- g^ 

d3q 



(27r) 



5* 95 



5* 



9v4 



g 



dA 1 



5 dB 1 dB 



and 



5*5* - 5^ 
d3q 



(27r)3 



d^q 

(27r)3 
5* dB 



g^g^ 



5* 



2 dpq 



2 

(1 + 2^^) , 



(l + 2n^) 



(4.32) 



dA 



5 



9A 



1 dA 



g^g-f - 5^ 5</'m 5*.9<[' - 5^ ^V'm 2 9/^4 



(1 + 2n^) 



5 



as 1 



(l + 2nB) 



(4.33) 



We now make use of the expressions for the condensate densities at finite temperature, Eqs. (|4.29() and (|4.30|l . with 
Eqs. H4.24II and 14.2511 together with the self-energies expressions Eqs. H4.26|l and (|4.27|l . to express Eqs. I|4.32|l and 
(|4.33l) completely in terms of the temperature dependent condensate densities (e.g., in terms of 0,„ and ipm) instead 
of the tree-level (T — 0) condensate densities 4>q and iJ^q. This process is analogous to the one used to obtain Eq. 
(|4.19l) for the one-field case. After some straightforward algebra, this then results in the coupled equations expressing 
Ipm and (pm in terms the total densities p^. and 



d^q 

(27r)3 



OA , „ , dB , „ ; 
(l + 2?i^) + - — (1 + 271^) 



dp<s. 



(4.34) 



and 
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Using Eqs. 



and (|3.21l with /i 



(27r)3 



OA , . dB , 

(l + 2n^) + - — (1 



fii, we obtain that 



(4.35) 



with 



dA 


1 


dijL^ 


"il 


dA 


1 




^AA 


dB 


1 




^4B 


dB 


1 







-G$) + 



1 

1 




— — f r 1 1 ii ,T* f T iT* ft ^ 1 r ^ 1 - 

V. ' i \ y/ ^ ' y/ J. J. 1^ Vw-* 1^ / 


- 2o2(i2 W;2 


4A 






^mV'mG.I.G* 


1 




+ G>i/) (^H^iG-qi — -ff$G<j)) 


+ 2520^V'^G$ 


AA 


V 


(iJvpGvi, ~i7<E,G$)%452, 


^m'0mG*G* 


1 




+ G$) (^H^G'ii — H^G^^ 


- 2520^^^G* 


45 


V 


{Hq,G^ - Hq.G^y +4:g^^ 




1 




+ G\[/) (-ffvi>Gvi> — -ff<i)G<i>) 


+ 2^20^ V^G$ 


45 


V 


^(i?^,G^, -i7$G$)V 4^2 


^m'/'mG<i,G^, 



G^ii 
G<i> 



2™^ 

2m$ 
2m$ 



(4.36) 



(4.37) 



The coupled equations H4.34|l and H4.35|l seem very messy. But still we can obtain a few analytical results from it 
and perform some qualitative discussions about the phase diagram {p^, p^,T). for the coupled system. For instance, 
if we consider an equal mass system, m$ = = m, the expressions for A and B simplify to 



A^ 
B^ 



q / q~ 

2m \ 2m 



2m V 2m 



(4.38) 



with 



9'S'IPr, 



± 



2 i2 rz 
9 Vm^Pr, 



1/2 



(4.39) 



Using Eqs. (|4.38|) and (|4.39|) in Eqs. (|4.34|l and H4.35|l we can compute the quantum depletion terms (the T = 
terms) appearing in the coupled system of equations, Eqs. (|4.34() and H4.35|l . Like in the one field case, we use again 
dimensional regularization to compute the zero temperature momentum integrals in Eqs. 1)4.34(1 and 1(4.35(1 to obtain 
the results. 



f d^<l 


/ dA 


^ dB\ 


{2ma+f/'^ + (2ma_)3/2 


J {2nY 




dp<i, ) 


127r2 



14 



(2ma+)3/2 - {2ma-) 



3/2 



2 2 
2 



1/2 



127r2 



V 87r2 
(—)■ 



5*— ( 



1/2 



- (2ma_)i/2 (4.40) 



and 







/ dA 


dB \ 


(2ma+)'V2 + (2ma_)3/2 




(27r)'i 






127r2 



(2ma+)3/2 - (2TOa_) 



3/2 



V 2 



g*? 



2 2 
2 



1/2 



127r2 



m 



V 8^2 
/ m 



(2ma+)^/2 + (2TOa_)^/2 



1/2 



(2ma+)^/2 - (2TOa_)^/^] (4.41) 



From cither Eq. H4.40|l or Eq. (|4.41|) . it can easily be checked that for g = we re-obtain the result (|4.2(J|) for either 
the ^ or the $ fields. Note also that by taking g = Eqs. H4.34|l and H4.35|l decouple and we recover the one-field 
expression H4.19|l for each of the fields individually. It is interesting to point out that, at this level of approximation 
we are considering, the above equations show that if any of the fields go above the transition point (either (j)m = or 
"fprn = 0) the two equations H4.34|l and H4.35|l decouple, becoming independent of each other, since the cross-coupling 
term in (|4.36(l always appears multiplying both (j>m and 4'm- As far as SNR/ISB are concerned one may conclude, 
based on the above equations, that these phenomena do not arise for this theory since the cross coupling always 
appears as g^ so the relevant physical quantities are insensitive to the sign of g. 

Equations H4.34|l and H4.35|l are the two-field analogous of the one-field case, Eq. (14.191) . Note that here it is more 
convenient to express the resulting expressions for the densities completely in terms of and tpm ■ In the one-field 
case both the condensate density and the effective chemical potential, for definition, vanish at the critical point. In 
the two-field case, as Eqs. ()4.29|l and H4.30|l show, at the critical points for <I> and 'J we have (f>m{T = T$) = and 
fpm{T = T\s,) = 0. However, the same does not necessarily (as it should not, actually) happens with the effective 
chemical potentials, /!$ and The critical temperature for transitions in the $ and \1/ directions are to be determined 
from the numerical solution of the coupled equations H4.34|l and H4.35|l . For this, it is useful to express Eqs. (|4.34|l 
and (|4.35|) into dimensionless quantities, like in Eq. (|4.22|l for the one-field case. We start by making the definitions 
pxs, — Op^, nixi, = TO$ —171,0^^— 7011(0$ and at — {ni/ piY^^ where ni {i — 5", $) are the diluteness parameters for 

and <i>. It then follows that gi — {At: / 'm){ni / piY^^ , g2 — 7(87r/m)2(rt,p/p,p)i/3(n$/p$)i/'^, whereas the boundness 
condition now reads just 7 < 1/4. As in the monoatomic case, = q^/{2mT). Other useful quantities are the 
dimensionless temperatures Ti = T/Ti^c, where Ti^c = (27r/TO)(/9i/C(3/2))2/3 represents the critical temperature for 

the monoatomic case, and the dimensionless integral measure, (l / pi)[d^q/ {2Tr)^] = {1/ C,{2)/2)){A/ ^jTT)Tf^'^x^dx. Then, 
Eqs. H4.34|l and f I4.35|) can be written as 



1 

2^ 



(2^ 



dA 
dp.i 



dB 

dp.i 



1 



dh 



(2^)3 \^p^ 



dA 



dB 

dm 



(4.42) 



where the temperature independent (quantum depletion) part is given by 



15 



d^q ( dA dB 



J {2tt)^ \dfj,^ dfj,^ 



+ 



- -(n^/W)[(2a+)'/' + (25-)^/'] 



1/3- 



1/3 1/3- 



1/2 



{ K^'p^.c - ni/^p$,e/e2/3]2 + 16(7/02/3)n^/3^i/3 -^^^ -^^^^ | 

X [(2a+)i/2 _ (2a_)V2] , (4.43) 



where 



a± = (47r) ^ n^^^p*,c + ^^^tI ± 



Then, using p$ = /0 one gets 



X 2 

1/3- l/3P*,c\ , i« 7 1/3 1/3- 

P*,c - ^273 ) + ■'^"^273'^* /^*.cP4>,c 



1/2^ 



(4.44) 





/ 


dB \ 


/ (27r)3 


\dp^ 


dp<s,) 



127r2 



[(2ci+)3/2 + (25_)3/2] 



i/'p*,c - ni/=^p$,e/e^/^ [(2a+)3/2 - (2a_)3/2] 



,1/3- 



^«/^p*,e)[(2a+)V2 + (2a_)V2] 

1/2 



+ 



{ K^'p^.e - ni/^'p^.e/^^/^l^ + 16(7/^2/3)n^/\J/3p^,,p$,,) } 
X [(2a+)^/2 - (2a_)i/2] . (4.45) 



Let us now use = T/T^^c as a reference temperature to define the temperature dependent parts, starting with 



p^l J {2-K 



d^q ( dA 



riA 



dB 
dp-a 



riB 



Jo 



dxi 



dA{x ) 
dp<s 



dB{x) 
dp<s 



(4.46) 



where the dimensionless Bose factors are ua^x) = 1/(6"^^^' — 1) and nB{x) = l/(e^(^' ~ 1); with 

1/2 



A{x) = X 



5-C(3/2)^/n 
27rTW J 



(4.47) 



and 



V 27rf* y. 



(4.48) 



Let us further define 
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.x2 + 2C(3/2)2/3nf 



(4.49) 



and 



Then, 



and 



Six) 



.2 + 4C(3/2)^M^^ 



/c(3/2)y/'^i/3P*^ 



H{x) = 



^2^fi4 C(3/2)'/' 1/3 1/3- ~ 
g{x) +647 — — — P^,cP'S>,c^ 



1/2 



^*(a;)g(.T) + I67 



C(3/2)^/^ 1/3 1/3. . 



(4.50) 



(4.51) 



(4.52) 



as 



1 



2B{x) 

Finally, let us write down 



J-^ir{x) 



1 , ^^^r C(3/2)^/^ 1/3 1/3- - x^^ 

I ^*(a;)g(a;) + I67 — ^ — P^f,cP<s>,cj:^ 



(4.53) 



1 f d^q f dA dB 



T',?/^6l / dx 



(dA{2 



using 



dB{x ) 
dpi 



(4.54) 



JF$(x) = 

and the above definitions for Q and Ti.. Then, 



a;2 + 2C(3/2)2/3ni,/^-^ 



(4.55) 



and 



a4 



dB 
p<s> 



2Aix) 



J-^(x) 



1 



nix) 



^'S>{x)g{x) — 167 



C(3/2)4/3 1/3 1/3, 



92/3 



2B(x) 



J^$(a;) 



7^(a;) 



•?^*(a;)^(a;) — I67 



C(3/2)4/3^^,/3 ,/3^ ^ ^2 



612/3 



'T2 



(4.56) 



(4.57) 



We have carried out a careful numerical analysis finding that the cross coupling (now characterized by the dimen- 
sionless parameter 7) does not affect the critical temperature that each gas observes when the system is uncoupled. 
As shown in Fig. |21 when the temperature is increased each type of condensate returns to the symmetric (gas) 
phase at a critical temperature whose value coincides with the one obtained in the one field case (see Eq. H4.21(l 'l. 
That is, the two distinct critical temperatures displayed in Fig. |31are insensitive to 7 (i.e. g) and seem to depend 
separately on each density ^ that are now related by 0. Although the critical temperatures for each gas in the coupled 



* Actually, it is very plausible that the critical temperatures should also depend on the different masses as in Ea. l4.2H . However, here 
we had to use the approximation ~ = m which does not allow us to fully confirm this fact. 
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FIG. 3: The quantities pc,* and pc,* plotted as functions of the reference dimensionless temperature for 7 = 1/6, 6 = 1.5 
and = n# = 10~^ The solid lines correspond to the coupled case whereas the dashed lines are the corresponding curves for 
the one field case. Note how the observed Tc's agree with the value predicted for the one field case, Eq. 14. 2H . The symmetry 
0(2) referring to the less dense gas, <I>, is restored at a lower temperature — 0.763 (which corresponds to = 9T<s, = 1) 
than the one observed by the other specie, ^I" (T* = 1). 

case coincide with the ones observed in the monoatomic case the temperature dependence of the quantities pd are 
influenced by the cross-coupling. 

To investigate how the cross coupling influences the general temperature behavior of pc,i, we consider a large, a small 
and the null value for 7. Figure ^ shows the situation for the ^ sector displaying the fact that, for < Tij. < 1, the 
density pc,^ assumes smaller values for larger values of the cross coupling. However, as noticed before, the critical value 
for the temperature corresponds to the monoatomic case that in fact is the critical temperature for a non interacting 
gas. This is expected since at this (one-loop) level of approximation the Hugenholtz-Pines theorem washes out all 
nontrivial contributions to Tc (see Ref. so that our results for this quantity become trivial. Nevertheless Fig. 

21 suggests that the critical temperature value may be influenced by the cross coupling interaction in a computation 
that includes higher corrections. 




FIG. 4: The quantity Pc,* as a function of the dimensionless temperature T\i> for 9 — 1.0 and nip = n$ = 10"^ The solid line 
corresponds to a strongly coupled case, 7 = 2/9, whereas the dotted line corresponds to a weakly coupled case, 7 = 1/10. The 
dashed line corresponds to the uncoupled case. 
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V. THE SYMMETRY RESTORED, HIGH TEMPERATURE PHASES 

Lets us now study the two-coupled field system in the symmetry restored phase. In this case the spectrum for both 
$ and ^ field changes, since (j)m and 4'm both vanish. The functions and G$ defined by Eqs. 14.12|l and 14.13|l . 
and similarly for and G^, become 



2to$ 



(5.1) 
(5.2) 



where = ^$ — and /iip = /iip — Siit,!^, with the self-energies S, in the normal phase and at one-loop order, 

given by 



^01,01 ^02,02 




exp 



2 



exp 



f3 



2m a 



(5.3) 



and similarly for 



/^^jnN 3/2 

S*,* = ^^i,->pi = S^2,V'2 = ( ^_ ) Li3/2 (e''''*^ 



V 27r 



2ti 



3/2 



Li3/2 (e'^^*) 



(5.4) 



where Li3y2(-2^) is a polylogarithmic function, 



(5.5) 



The pressure in this case, P(T, /i^,) — V^ff (T, 00 = 0, ■00 = 0), becomes just 



P(r, ^*) 



1 (m^T\ 



3/2 



Li 



13/2 



(5.6) 



where the two last terms in Eq. H5.6|l come from the sum of two-bubble vacuum diagrams made of the quartic 
self-interaction and cross-interaction vertices with the two additional terms, proportional to the self-energies, added 
to the original Lagrangian and regarded as additional interaction terms (the last term in Eq. (|4.15|) and the similar 
contribution for the components of ^P). Using Eqs. (|5.3|) and (|5.4|) in H5.6|l . we obtain 



P(T,^$,/i*) = 



+ 9. (^)'Li3/2 (e'^^*) +5* (^)'Li3/2 (e'^^*) 
+ ff(m*m*)3/2 (^^^'03/2 (6/3^^)^3/2 (e''^*) 



(5.7) 
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From Eq. H5.7() we can now fix the chemical potentials from the expressions giving the densities, 



and 



where 



dP fm^T 



3/2 



Li 



3/2 



(3fj.il 



+ 



mqjT 
2n 



3/2 



Li3/2 (e'^'^*) 



(1 + 2.g^A^)- h gA^ 



dP /m^T 



V 27r 



3/2 



Li 



3/2 \e 



(1 + 



Li 



3/2 (e^'^-) 



;i + 2g$A<i,) 



dfH 



gA^ 



A, 



3/2 



Lii/2 (e'^^O 



The derivatives involving the chemical potentials in H5.8() and H5.9() , are defined by 



(5.8) 



(5.9) 



(5.10) 



9^$ 
9/2$ 



1 O /I /I 

1 — 2g$A$— ff^*" 



. 9p,p 
-2(7,p^>i,— (7^< 



1 — 295,^15, 
-2g*A$ 



9/i$ ' 



9/2 <!, 



5^* 



9/i* 



(5.11) 



S/iif 9p,i, 

Eq. H5.11|l represents a set of equations for the derivatives of the effective chemical potentials that can be easily solved 
and then the results substituted back in Eqs. (|5.8() and H5.9|l . The resulting expressions are just 



2n 



3/2 



Li3/2 (e'''^*) 



and 



P* - ( ^ I 1^13/2 



Li.v. (e''^*) 



(5.12) 



(5.13) 



Eqs. (|5.12l) and H5.13|l can also be used in the equations (at the one-loop level) defining the self-energies in the normal 
phase, Eqs. H5.3|l and (|5.4(l . Then, by expanding Eqs. (|5.12() and (|5.13(l in the high temperature limit {p,i/T <^ 1) one 
obtains the expression for the chemical potentials in terms of the densities (after inverting Eqs. (|5.12|) and HS.l^fl l. 



/i$ ~ 2g$p$ + gpq, 



P*-C(3/2) 



(5.14) 



and 



/i* ~ 2g^p^s, + gp<s, 



T 
An 



( 27r 



3/2 



P* - C(3/2) 



(5.15) 



The above results, for the uncoupled two Bose gas (g = 0) can be shown to agree with the results of Ref. obtained 
for the one Bose gas case and similar analysis that also follow form our results (|5.14|) and H5.15|l . 
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VI. CONCLUSIONS 

We have considered a nonrelativistic model suitable to describe a system of homogeneous dilute Bose gas composed 
by two different types of atoms. A survey of the literature [T3 | shows that there is a growing interest in this type of 
systems. In general, coupled systems show a richer phase structure in comparison with uncoupled ones due to the 
presence of a cross coupling 0,l2jlSE|- For example, coupled nonrelativistic systems under the influence of external 
fields, represented by one body terms, may exhibit re-entrant phases 0i 0- The appearance of such phenomena 
depends on the sign of the inter-species coupling, being completely ruled out in single species models. Apart from this 
fact, one also observes that the uncoupled and coupled models produce different values for critical quantities such as 

In this work, our main motivation was to check how the presence of a inter species coupling would affect the 
qualitative and quantitative behavior of the transition regarding homogeneous coupled Bose gases. With this purpose 
we have evaluated the effective potential at finite temperature in a nonperturbative fashion to one-loop. Due to the 
complexity of zero temperature contributions, the complete evaluation was only possible in the approximation where 
both atomic masses are approximately the same. Our results show a dramatic difference concerning the case studied 
in Ref. where the one body term represents, e.g., external fields and the Bose gas case, considered here, where 
these terms represent chemical potentials. Our present results exclude the possibility of exotic transition patterns such 
as inverse symmetry breaking and re-entrant phases arising in models relevant for BEG. This is rather satisfactory 
since, intuitively, one expects that the BEG transition for coupled gases should also observe the same simple pattern 
observed by monoatomic gases. That is, the system smoothly goes from an unsymmetric phase to a symmetric phase 
as soon as the critical temperature is reached. This is so because the value of cross coupling always appears squared 
eliminating the possible occurrence of ISB. Then, our second concern was to check the numerical values for the critical 
temperatures at which symmetry restoration occurs. Numerically, we found that the Tc values are insensitive to the 
coupling and coincide with the monoatomic results (which at this level of approximation corresponds to the standard 
ideal gas result). This can be understood in the view of the Hugenholtz-Pines theorem, in which case our one-loop 
nonperturbative approximations cannot probe the effects of interactions at T = T^. Nevertheless, in the region 
< T < Tc (where our approximation is more reliable) the curve describing the transition is seen to be influenced 
by the cross coupling. Taking the two species as having different densities, we have started at T = increasing the 
temperature and observing a first transition in which the more dilute specie always reaches the non-condensed phase. 
At this first critical point the system decouples, the denser specie remains in the condensed phase while the less dense 
acts as a thermal cloud. Then, a second transition occurs with the denser specie reaching the gas phase exactly at 
the Tc observed in the corresponding monoatomic case. If the two species have the same density the unbroken phase 
is reached at once, which is not surprising, remembering that we took the atomic masses as being approximately the 
same. In summary, as in the monoatomic case and due to the Hugenholtz-Pines theorem, only the shapes of the curves 
describing the thermal behavior of the effective potential minima are sensitive to the numerical values of the couplings. 
The cross coupling, in particular, has a non negligible effect in this parameterization hinting in the possibility that a 
shift will appear if one computes, at higher orders, the critical temperatures for interacting uncoupled and coupled 
Bose gases. 

Though in this work we have considered a rather idealized case, in comparison to real experiments performed in BEG, 
by only working with a homogeneous model and neglecting any non-homogeneity effects, we hope that the approach 
developed here may still be useful in the analysis of the thermodynamics in realistic Bose-Einstein condensation 
experiments with coupled atomic gases. The results and field theory methods used here could be considered, for 
instance, to be applicable to trapped atomic gases in the central region of wide traps. 

Finally we remark that the coupled field model studied here, Eq. H2.1|l . could actually be viewed as the nonrelativistic 
limit of a two complex scalar fields $ and 5*, both with conserved charges, with chemical potentials i = In 
the nonrelativistic limit, the chemical potentials appearing in Eq. H2.1|l should then be identified with /i^ = /i^.i — 
which is the correct identification of a nonrelativistic chemical potential ^] . Some of the details of the nonrelativistic 
limit of the corresponding relativistic action of a two-scalar field model were already given in Refs. a,|ci, while for 
the one field model this was previously discussed in |22| | . The calculations and analysis performed in the present work 
could then be easily extended to the relativistic problem and the combined effects of temperature and finite densities 
for the corresponding phase diagram of the model be studied. This would be of particular relevance for studies related 
to the early universe phase transition problems as well as the current heavy-ion collision experiments. We will report 
on the extention of the results presented in this work to the relativistic model in a forthcoming publication. 
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